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The NPDGamma y-ray detector has been 
built to measure, with high accuracy, the 
size of the small parity- violating asymme- 
try in the angular distribution of gamma 
rays from the capture of polarized cold 
neutrons by protons. The high cold neu- 
tron flux at the Los Alamos Neutron 
Scattering Center (LANSCE) spallation 
neutron source and control of systematic 
errors require the use of current mode 
detection with vacuum photodiodes and 
low-noise solid-state preamplifiers. We 
show that the detector array operates at 
counting statistics and that the asymme- 
tries due to B4C and Al are zero to with- 
in 2 x 10~* and 7 X 10~^, respectively. 
Boron and aluminum are used throughout 
the experiment. The results presented here 
are preliminary. 
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1. Introduction 

The NPDGamma experiment is under commission- 
ing at the Los Alamos Neutron Scattering Center 
(LANSCE). It is the first experiment that was designed 
for the new pulsed high flux cold beam line, flight path 
12, at LANSCE. NPDGamma will determine the very 
small weak pion-nucleon coupling constant f^ in the 
nucleon-nucleon interaction [1-5]. This coupling con- 
stant is directly proportional to the parity- violating up- 
down asymmetry Ay in the angular distribution of 2.2 
MeV gamma rays with respect to the neutron spin 
direction in the reaction n + p ^ d + y. The asymme- 
try has a predicted size of 5 x 10"^ and our goal is to 
measure it to 10 %. 

The small size of the asymmetry and the high pro- 
posed measurement precision impose heavy require- 
ments on the performance of the beam line and appara- 
tus. It is necessary to achieve high counting statistics 
while at the same time suppressing any systematic 
errors below the statistical limit. The detector array was 
designed to satisfy these requirements [6,7]. 

The experiment uses an intense, pulsed (20 Hz), cold 
neutron beam at LANSCE [8]. The beam is transverse- 
ly polarized by transmission through a polarized ■'He 
cell. Three ^He ion chambers are used to monitor beam 
intensity and measure beam polarization through trans- 
mission ratios. A radio frequency spin flipper is used to 
reverse the neutron spin direction on a pulse by pulse 
basis. For the production experiment the neutrons will 
capture in a 21 L liquid para-hydrogen target. The 
gamma rays from the capture are detected by an array 
of 48 CsI(Tl) detectors. The entire apparatus is located 
in a homogeneous 1 mT magnetic field to maintain the 
neutron spin downstream of the polarizer A detailed 
introduction to the experiment is given in [1]. 

During commissioning the radiative neutron capture 
on various target materials was investigated to look for 
any false asymmetry which may enter as a systematic 
effect while taking data with the hydrogen target. The 
measurements concentrated on materials that can be 
found within the experimental apparatus and which are 
interacting with the neutrons. The targets used included 
Al, Cu, In, ^Li, and B4C. 

Boron is used throughout the experiment, for neutron 
shielding and to collimate the beam. Aluminum is used 
in most of the equipment and the beam encounters sev- 
eral millimeter of it, primarily in the windows of the 
hydrogen target. It is therefore necessary to establish 
the size of the y-asymmetry due to neutron capture on 
both of these elements. A preliminary result on the alu- 
minum asymmetry is presented in [1]. The preliminary 



asymmetry result for the B4C data taken so far is pre- 
sented here. Further analysis on all targets is currently 
underway. 

To establish that the array operates at counting statis- 
tics, it is desirable to conduct a neutron capture experi- 
ment in which a single neutron capture produces at 
most one outgoing gamma ray with a known energy. 
The reaction '''B(n,a)^Li produces an excited ^Li nucle- 
us with a branching ratio of 94 %. The excited nucleus 
decays, emitting a 0.48 MeV gamma ray, 100 % of the 
time. As a result, the width in the data taken with a B4C 
target can be compared to the width expected from neu- 
tron counting statistics. 

Based on LANSCE FP12 beam brightness measure- 
ments and Monte Carlo calculations [8], the predicted 
average gamma rate into a single detector is 20 MHz. 
At that rate, pulse counting is not possible, given the 
long decay time of the scintillation light pulses in Csl. 
The detectors are therefore operated in current mode. 
The scintillation light from the Csl detectors is convert- 
ed to current signals using vacuum photodiodes (VPD), 
and the photocurrents are converted to voltages and 
amplified by solid-state electronics [7]. 

The 48 detectors are grouped in rings of 12 detectors 
each (Fig. 1) and arranged in a cylindrical pattern 
around the neutron capture target used. The detector 
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Fig. 1. A ring of detectors and one up-down pair, as seen in beam 
direction. B is the magnetic holding field defining the direction of 
the neutron polarization. 
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signals are sampled by 1 6 bit analog-to-digital convert- 
ers (ADCs) and integrated over a 0.4 ms period, defin- 
ing the time bin width at which the detector signals are 
stored. For each neutron pulse, a total of 100 time bins 
per detector is stored, followed by a 10 ms break before 
the next frame of neutrons arrives. This 50 ms period 
defines a macro pulse of data. 

The remainder of this paper gives a short description 
of the counting statistics and asymmetry tests per- 
formed and the relevant preliminary results obtained. 



2. Achieving Counting Statistics in 
Current Mode Operation 

In current mode detection, counting statistics appears 
as the RMS width in the sample distribution, due to the 
fluctuation in the number of electrons produced at the 
photocathode of the VPD, as a result of the quasi 
instantaneous amount of energy deposited in the Csl 
crystal. 

During beam on measurements, the shot noise RMS 
width is given by [9] 



^4ot = V2^VZ"' 



(1) 



where q is the amount of charge created by the photo- 
cathode per detected gamma-ray, / is the average pho- 
tocurrent per detector andf^ is the sampling bandwidth, 
set by the 0.4 ms time bin width. 

The boron target used in this measurement consisted 
of a 1 cm thick 1 5 cm by 1 5 cm sheet of sintered B4C 
located in the center of the detector array. The B4C tar- 
get is black for cold neutrons. The beam was collimat- 
ed down to a 10 cm diameter cross-section, so that the 
target completely covered the beam. Three consecutive 
measurements were made. Three 500 s long pedestal 
runs, without beam, one 500 s data run with beam inci- 
dent on the target and one 500 s background run with 
beam but without the target installed. 

Because of the single gamma released in the 
'''B(n,a)'Li capture reaction and Eq. (1) a direct com- 
parison can be made between the RMS width expected 
from counting statistics and the RMS width seen in the 
data stream. Several factors contribute to the overall 
width seen in the data stream. Contributions from elec- 
tronic noise and signals from activation decay in the 
various materials are measured during beam-off, 
pedestal runs. Contributions from radiative capture and 
scattering on materials other than the intended target 
are measured by performing background runs. Both of 
these contributions can be removed by subtracting their 



RMS widths in quadrature from the overall width. In 
addition, any fluctuations in the magnitude of the sig- 
nal over the course of the measurement will also con- 
tribute to the width. To remove those contributions, the 
data and background runs were normalized to the beam 
current and separate variances were calculated for each 
group of 20 macro pulses and then averaged. For all of 
these calculations a single time bin, corresponding to 
the peak signal (£" = 14 meV), was used in each macro 
pulse. 

Figure 2 shows the RMS width for a typical detector, 
as seen at the preamplifier output. For this detector, 
using Eq. (1), the measured values for q and / and the 
gain of the preamplifier circuit, the RMS width expect- 
ed at the preamplifier output is 5.7 mV ± 0.3 mV. The 
error on the expected width is dominated by the accu- 
racy to which we know the efficiency (number of pho- 
toelectrons per MeV) of the detector Additional contri- 
butions to the width, that are not currently accounted 
for, include fluctuations in the amount of energy a 
gamma-ray actually deposits in the detector (currently, 
it is assumed that each photon deposits its full energy) 
and background due to neutrons that are scattered off 
the target and subsequently capture in other surround- 
ing materials. An effort to carefully model these effects, 
using Monte-Carlo methods is currently underway. 
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Fig. 2. Counting statistics analysis results for a typical detector. The 
RMS width due to counting statistics is compared to the width seen 
from pedestal runs (electronic noise). The fit to the beam on, target 
in data histogram shows an RMS width of 6. 1 mV ± 0.04 mV. 



Using the known proton current and an appropriate 
Monte-Carlo model to simulate the transmission of the 
neutrons through the guide and the experiment, the 
number of neutrons per pulse expected to capture on 
the B4C target, for the peak energy, is about 8 x 10'' 
ms~\ Considering the gamma attenuation in the target 
as well as the solid angle for the detector, the number of 
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photons into a detector, per time bin is 3 x lO'*, so that 
the relative width is 77^= 0.6 %. Using the correspon- 
ding vohage, one expects to see out of the preamphfi- 
er, the RMS width expected due to neutron counting 
statistics is =5 mV ± 1 mV. The error is due to the accu- 
racy to which we know the detector efficiency (=5 %), 
the estimate on the accuracy in the number of neutrons 
capturing in the target (=5 %), but mostly due to an 
error estimate on the number of gamma rays that will 
deposit all of their energy in the detector (=20 %). 



3. The Parity Violating y-asymmetry in 
"B(n, y)"B and "Al(n, -^'^Al 

Due to the 1 mT holding field surrounding the exper- 
imental apparatus, the neutrons are polarized vertically 
after leaving the ■'He spin filter While taking hydrogen 
data, the parity violating asymmetry in n-p capture is 
therefore seen in a difference of the number of y rays 
going up and down. The cross section is proportional to 
1 +A^ cosO, where is the angle between the neutron 
polarization and the momentum of the emitted photon. 
Accordingly, when investigating possible false asym- 
metries, i.e., those not due to n-p capture in the hydro- 
gen target, we look for them in the up-down direction 
with respect to the neutron polarization axis. 

The geometry of the B4C target was described earli- 
er The aluminum target consisted of 30 sheets of alu- 
minum, each of which was =1.6 mm thick. Each sheet 
was separated by 5 mm to reduce gamma attenuation in 
the target. The total target thickness was 5 cm and its 
cross-section was 9 cm by 9 cm. The target was locat- 
ed in the center of the array. 

During the first phase of commissioning, approxi- 
mately 1.5 hours of B4C data were taken and the results 
show a zero asymmetry at the 2 x 10"^ level of accura- 
cy. The aluminum asymmetry was measured to be zero 
at the 0.73 X 10"^ level using about 40 h of data. 

In calculating the asymmetry for the entire period, 
one asymmetry was calculated for each detector pair 
and time bin and over any valid sequence of eight 
macro pulses with the correct neutron spin state pattern. 
A so-called valid eight step sequence of spin states is 
defined as (TsLiTiTTi). Using this pattern suppresses 
first and second order gain drifts within the sequence. 
Some eight step sequences may become invalid due to 
missing pulses in the DAQ. 

A pair of detectors is defined as shown in Fig. 1 . If 
we let (f/t, Ui or D^, Di) denote the sum of all four sig- 
nals with the corresponding spin states in a given spin 



sequence, then the detector pair asymmetry for a given 
time bin is calculated from 



^.w(0 = 






(2) 



This method of calculating the asymmetry allows for 
automatic subtraction of pedestals and any additive 
background. Signal fluctuations that are not correlated 
with the switching of the neutron polarization direction 
will average out and don't contribute to any asymme- 
try. Possible false asymmetries due to electronic pickup 
and possible magnetic field induced gain changes in the 
detector VPDs have previously been measured and are 
consistent with zero to within 5x10"' [6]. 

The asymmetry was then divided by the neutron 
polarization and the cosine of the corresponding detec- 
tor pair angle. For each detector pair the resulting time 
bin physics asymmetries within a sequence were com- 
bined by calculating an error weighted average over the 
time of flight selected for analysis (13 ms through 34 
ms). A 48 % ■'He polarization was calculated from the 
monitor ratios and was seen to be stable over the ana- 
lyzed runs. 

The spin-flipper efficiency was taken to be 100 % 
over the analyzed time of flight range. The resulting 
physics asymmetry for each valid eight step sequence 
was then histogrammed to produce the total asymme- 
try. The results for the asymmetry are shown in Fig. 3 . 

The asymmetry on aluminum was calculated using 
the same method. The results are shown in Fig. 4. 
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Fig. 3. Up-Down gamma-ray asymmetry for -1.5 h of B data. 
From the fit, the asymmetry is (-0.2 ± 1.9) x 10"^. 
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Fig. 4. Up-Down gamma-ray asymmetry for -40 h of Al data. 
From the fit, the asymmetry is (-0.95 ± 0.73) x 10"*. 



4. Conclusion 

The NPDGamma experiment successfully complet- 
ed a commissioning run in April 2004. The Csl detec- 
tor array was demonstrated to be sufficiently low noise 
to allow operation at counting statistics. The array was 
used to perform measurements of asymmetries in neu- 
tron capture in gamma-ray emission following capture 
of '"B and ^'Al. Given the thickness of the targets used 
and the amount of actual material that is estimated to 
interact with the neutrons, about 10 % of the signal 
seen with the hydrogen target will be due to aluminum 
while about 1 % comes from boron. The amount of alu- 
minum data analyzed corresponds to about 1/3 of the 
data taken during the commissioning run. Taking these 
factors into account and using all collected data, we 
will be able to show that neither material is going to 
contribute a false asymmetry in the background signal 
of the NPDGamma apparatus. The Csl detector array is 
fully commissioned and prepared for production data 
taking with the hydrogen target. 
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